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Abstract 

We discuss the approach of effective field theory on a <i-dimensional Euclidean space in a scalar 
theory with two different mass scales in the presence of flat surfaces. Then considering Dirichlet 
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1 Introduction 



The quantum field theory of self-interacting scalar fields has long served as a laboratory for de- 
veloping methods of analysis that can be applied to theories of more direct physical interest. Our 
purpose is to investigate the Euclidean field theory of two interacting scalar fields where after the 
construction of an effective theory for the light modes, in some limit, there is a decoupling between 
the light and heavy modes as stated by the Appelquist-Carrazzone theorem M. Aditionally we 
impose boundary conditions on the resulting theory in order to study finite-size effects and the 
renormalization program in systems where the translational symmetry is broken. The interest 
of the literature in the study of quantum fields in the presence of boundaries appears after the 
problem investigated by Casimir more than fifty years ago @|. A complete review of this effect 
can be found in refs. @ @1 @- 

In 1948, Casimir showed that neutral perfectly conducting parallel plates in vacuum attract 
each other, the effect can be interpreted as follows: the distortion of the vacuum fluctuations 
of the quantized electromagnetic field due to the presence of metallic plates renders the zero- 
point energy of the field into a measurable quantity. In the absence of any classical background, 
the renormalized vacuum expectation value of the Hamiltonian operator can be correctly defined 
by the Wick-ordered product. The main support for this procedure comes from the proof that 
for a relativistic field theory, the vacuum expectation value of the stress-energy tensor should 
vanish in order to ensure that the realization of the Poincare generators in terms of the fields 
of the theory satisfy the correct commutation relations 0. On the other hand, where external 
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fields or macroscopic structures are present, more elaborated methods must be used to find the 
renormalized vacuum energy of the quantized field while avoiding undesirable divergences. A 
procedure that identifies the divergent contributions to the vacuum energy, such as a cut-off or 
any analytic regularization method, followed by a renormalization is mandatory. The fundamental 
idea of the Casimir renormalization procedure is that although formally divergent, the difference 
between the zero-point energies of two different physical configurations can be shown to be finite. 

In quantum electrodynamics there is a standard argument used to support the implementation 
of a regularization procedure followed by renormalization to obtain the renormalized vacuum 
expectation value of the stress-energy tensor associated with the Maxwell field in the presence 
of boundaries. At high frequencies no real material is a perfect conductor. A wavelength cut- 
off corresponding to the finite plasma frequency must be included in the model. High energy 
modes are insensitive to the boundaries and only the low energy modes are affected by them. 
Consequently, in the study of quantum electrodynamics in the presence of conducting boundaries, 
in the generating functional for the n-point correlation functions this latter condition means that 
one may integrate out all the Fourier modes associated with the Dirac field and obtain an effective 
theory for the Maxwell field. 

We would like to implement the above discussed situation in a very simple model of two 
Euclidean self-interacting scalar fields. We will analyse a theory with two massive scalar fields 
with different mass scales and we will be interested in obtaining the effective action for the light 
field. One could theoretically envision a theory with two massive fields with different mass scales on 
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which the mass of the heavy field is smaller than the natural cut-off of the boundary, nevertheless 
in this paper we are not interested to discuss this situation. It is important to keep in mind 
that in order to construct an effective action that gives a correct description of the physics of the 
light modes in the presence of the boundaries, the Fourier modes associated with both fields with 
wavelength smaller than A" 1 for some cut-off A must be integrated out, restricting the space of 
functions that we are integrating over (note that we are always assuming a sharp cut-off). Since 
dealing with functional integrals with cut-offs in general models is quite complicated, we will limit 
ourselves to a heavy field with Gaussian functional integrals. We are studying finite size effects 
for the light modes in two steps. In the first, we integrate over the modes of the heavy field, 
obtaining an effective action for the light field. In the second step we are taking the limit in which 
the decoupling theorem is valid {m-i — > oo), to regard the effective action as the fundamental 
action and we are assuming boundary conditions over the remaining light field. A combination 
of different analytic regularization procedures and a renormalization are able to eliminate the the 
usual bulk and also the additional surface divergences that appear in the theory. The final result 
of our procedure is that we have the effect of the compactification of one dimension, breaking 
the full translational invariance of the original theory. In this situation we have an effective field 
theory of the light modes, with finite size effects. Of course, the region outside the boundaries is 
the union of two simple connected domains and the renormalization of the interacting field theory 
in such region must be carried out along the same lines that we have to use in the interior region. 
For simplicity we are only focussing our analysis to the interior region. 
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There are many papers in the literature discussing quantum field theory in the presence of 
boundaries or macroscopic structures. The calculation of the radiative corrections to the renor- 
malized energy density associated with the Maxwell field, in the presence of perfectly conducting 
plates, assuming no boundary conditions for the Dirac field was performed by Bordag et al J7|. 
Temperature corrections for this model were analyzed by Robaschik et al ||. Bordag et al also 
studied the leading radiative correction to the renormalized energy, assuming that the paralel 
plates are represented by delta function potentials ||. Using the approach of effective field theory, 
the radiative correction to the Casimir effect integrating out the fermionic degrees of freedom was 
examined by Kong and Ravndal and Ravndal and Thomassen [IIJ. A different approach was used 
by Falomir et al. JTTJ] . These authors studying scalar fields in the presence of a spherical shell 
used a sharp cut-off, assuming that the boundary is transparent for the heavy modes while for the 



soft modes they assumed Dirichlet boundary conditions. Actor |12| studied two interacting scalar 
fields in the presence of macroscopic boundaries assuming that only one of the fields satisfies clas- 



sical boundary conditions. Using the generalized zeta function method [ 13 1 the one- loop effective 



action was presented. More recently, Melnikov [|14]] investigated the low-energy effective action 
in a model with two scalar fields and also in quantum electrodynamics. Our treatment is very 
similar to the treatment developed by Melnikov in the study of the low-energy effective action for 
a theory with two interacting scalar fields. 

It is important to point out that the combination of effective field theory and finite size ef- 
fects can produce unexpected new phenomena. A well-known example of this situation is the 
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Scharnhorst effect f[q| . Studying quantum electrodynamics between perfectly conducting plates, 
Scharnhorst concluded that the speed of light normal to the plates, exceeds light speed on vacuum, 
while parallel to the plates light travels with its vacuum speed, i.e., with that of light propagating 
in unbounded space. Further calculations by Barton \W\ and also by Barton and Scharnhorst fT7|j 
confirmed the original result. The Scharnhorst effect is probably a consequence of the combined 
use of the effective Lagrangian in quantum electrodynamics and the presence of the plates, i.e., 
the Euler-Heisenberg Lagrangian density [18], derived many years ago and also rederived later 



by Schwinger [Tj|. As we already pointed out, the central issue of this effective Lagrangian is a 
derivative expansion of the photon effective action obtained by integrating out the fermionic field 
in the Maxwell-Dirac action. For an interesting discussion concerning the velocity of propagation 
of signals in different field theory models, see for example, JHJ . 

Finite size effects that do not break translational invariance in quantum field theory also have 
been extensively studied in the literature |H|] [fffi] . For translationally invariant systems, we 
can change from coordinate space to momentum space representation, the latter being a more 
convenient framework to analyze the divergences of the n-point Schwinger functions on which 
translational invariance is realized through conservation conditions. For systems where the trans- 
lational invariance has been partially broken (so there is still translational invariance along cer- 
tain directions) a more convenient representation for the n-point Schwinger functions is a mixed 
momentum-coordinate representation. Important references discussing the renormalization pro- 



gram in the presence of boundaries are the Symanzik [23] and also Diehl and Dietrich [24] papers. 
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In this paper we are studying the renormalization program in the presence of surfaces where a 
scalar field satisfies boundary conditions. We are interested in investigating a very simple model 
where we can construct an effective field theory for the light field on which the decoupling theorem 
holds. Further we consider an interacting Euclidean field theory of the light field in the presence 
of boundaries. We will consider a Casimir-like configuration where one of the coordinates, z, 
lies in the interval [0, L] imposing Dirichlet-Dirichlet boundary conditions and for the sake of 
completeness we will also enquire about Neumann- Neumann boundary conditions. 

The organization of the paper is as follows: In section II we introduce a simple model of two 
Euclidean interacting scalar fields. By integrating out the heavy modes associated with one of the 
fields, we are able to build the effective action for the remaining field. In section III we discuss 
a scalar field theory with boundary conditions due to the finite size of one of the coordinates 
and build the free two-point and four-point functions, both for Dirichlet-Dirichlet and Neumann- 
Neumann boundary conditions. In section IV we discuss the surface divergences of the one-loop 
two-point and also four-point function function. In section V we discuss the global approach, 
used to define the Casimir energy associated with a field in the presence of boundaries with well 
defined geometric shape.. Finally, section VI contains our conclusions. Throughout this paper we 
use h — c — 1. 
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2 The Euclidean functional integral and the effective ac- 
tion 



The goal of this section is to present a very simple model of two Euclidean self-interacting scalar 



fields where after the construction of an effective action p5| |26| |27| [ 28fl , and further imposition 
of the infinite mass limit for the heavy field, the decoupling theorem holds. For this purpose 
we start from a model with two different mass scales. We consider two real interacting massive 
scalar fields <fi(x) and <f2{x) with masses m 1 and m 2 , respectively and regard the field <f2{x) 
as the heavy field (m 2 >> mi). For this theory, we will rederive the Euclidean version for the 
Appelquist-Carrazonne theorem stating that for certain choice for the self-interacting (fi(x) part 
and in the limit m 2 — > oo, there is a decoupling in the effective theory. The only effects of the 
heavy field <£" 2 (a;) being a modification of the value of the renormalized mass and the coupling 
constant of the light field (p\(x). 

We start from the generating functional for the n-point Schwinger functions associated with 
two massive real fields in a (^-dimensional Euclidean space given by 

Z[j uj2 ] = M j [#l][# 2 ] e-SlVuwMsourceterms)^ (1) 

where [d<pi] [d<p 2 ] = I\ X €R d dipi(x)dip 2 (x) is an appropriate measure, S[tpi, y5 2 ] is the classical action 
associated with the Euclidean fields, and in the generating functional, M is a normalization. As 
usual, the n-point Schwinger functions of the theory can be obtained by functional differentiation 
with respect to the external sources ji{x) and j 2 (a;). Since in this section our interest is to construct 
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the effective theory for the light field, the introduction of the external sources in the functional 
integral is not important for our discussion. We consider the theory described by the following 
Euclidean Lagrangian density with two real scalar fields 

C(ifi, <p 2 ) = £o(y?i, ¥2) + Cint(vi, ¥2), (2) 

where the free part of the Euclidean Lagrangian density is given by 

1 11 1 

Cq[<pi, ¥2) = ^(^i) 2 + 2 m ?^? + 2 (^ 2 ) 2 + 2 m ^2, (3) 

and the interacting part is given by 

Antfal, ¥2) = Vfa) + y (¥l¥>2) 2 - (4) 

It is important to remark that the precise form of V(ipi) is not important for the construction of 
the effective action, but as we will see later in this section, the form of the V(ipi) is important to 
implement the decoupling theorem. 
The action of the model is given by 

5^1,^2] = J d d xC{<£i{x),(p 2 {x)), (5) 

and using Eq.@ and Eq.(f|) can be conveniently split up as 

S[<pi,(p2] = S[(pi{x)} + S V2 [tpi{x),<p 2 {x)}, (6) 

S , [</9i(x)] being the (/^(^-independent part of it. In order to obtain a derivative expansion of the 
effective action r e j/[<£>i] we have to assume 1712 » mi. As usual, the operators (—A + mf) -1 and 
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(— A + m^) -1 must be used to define the free two-point Schwinger functions of the fields <pi(x) and 
(fi2{x) respectively, and we must recall that A stands for the Laplacian operator in lZ d . According 
to the above, the free two-point Schwinger functions of both fields can be represented by 

G(x-y',rm) =<x\{-A + m 2 i )- 1 \y>, i = 1,2 (7) 

and they obviously satisfy 

(-A + m*)G(x - y; m % ) = 5 d (x - y). (8) 

To obtain an effective action for the light modes of the theory, we integrate out the heavy field 
<f2(x) in the functional integral and define the effective action of the light modes r e //[<^i] by 

e -r e// [ Vl ] = Ad^e-^vi.^]. (9) 

Using Eq.@ it is possible to write Eq.flj^) as 

e -r eff [ V i] = e -8[<pi] J ^ e -^ 2 [™]_ ( 1Q j 

The first step in our calculation is straightforward, since we have 

S ip2 [i fl ,if 2 ] = l^x^-A + m^ + ^O^) 2 ), (11) 

and using Eq. ([TT|) the functional integral appearing in Eq. (]10|) can be performed by means of 
Gaussian integrations, yielding 

e -r e/ /M = e-sirifato)-*, (12) 
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where we have 

0(x, y; m 2 ) =< x\0\y >= (-A x + vn\ + \^\{x))b d {x - y). (13) 

Consequently, the effective action for the light field ifi(x) is given by 

r e f f [<p 1 ] = S[<p 1 ] + ±trlnO. (14) 

Dropping a term that contributes trivially to the effective action T e ff[(px], we get 

r eff [tpi] = %i] + X -tr ln{l + X 2 (-A x + m*)-Vi). (15) 

There are many ways to evaluate the Fredholm determinant, defined by the above equation. Using 
a series expansion, it is possible to rewrite Eq.(^) as 



T eff fa] = Sfa) + - £ ^tr(A 2 (-A, + m 2 TVi)) fc , (16) 
or in a more compact notation 

oo 

k=i 

where each term of the series r( fc )[<^i] is given by Eq . (|16"D . Let us study the first non trivial 

contribution of this series, namely, the term k = 1 which in fact corresponds to a one-loop 
diagram. It is explicitly given by 

T (1) bi] = y J d d xG{x-x;m 2 )ip\(x). (18) 
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Using the Fourier representation for the two-point Schwinger function associated with the heavy 
field <p2{x), and defining a new coupling constant a = X 2 fi A ~ d , where \i is the usual dimensional 
parameter that appears in the dimensional regularization procedure, we readily obtain 

= WTT^m - i){m 2 ) d - 2 [ d d x ip\{x). (19) 



2(2^) d v 2' 

Now, the Gamma function T(z) is a meromorphic function of the complex variable z with simple 
poles at the points z — 0, — 1, —2, .. In the neighborhood of any of its poles z = —n, for n — 0, 1, 2, .. 
r(z) has a representation given by 

r(^) = ( " 1)w ] , +n(* + n), (20) 

n! + n) 

where Q(z + n) stands for the regular part of the analytic extension of T(z). Note that for odd 
dimensions r^ 1 )^!] is completely regular while in even dimensions (d = 2£ £ = 1,2, ...) there are 
singularities in the dimensional regularized quantity. By using the standard dimensional regular- 
ization prescription d = 2£ — e and since the Eq.([T9|) is quadratic in the field <fi{x), the divergence 
can be absorbed in the renormalized tpi(x) mass. Consequently, we define the renormalized mass 
for the light field as 



2 2 , am 2 



d-2 



m lR = m l + lA x d 
(47r) 2 



'-V)%- x 2 e x 
tf^T)\~e +n{ 2< 



(21) 



For the odd dimensional case there are no poles, but we have the same situation. We have thus 
shown that at the one-loop approximation, the first correction to the effective action given by 
r^ 1 ^^] that we obtain integrating out the heavy field <p2{x) is just a modification of the value of 
the renormalized mass associated with the light field. 
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We will now show that the second correction to the effective action given by r( 2 )[</?i] only 
modifies the value of the coupling constant of the field <fi(x). To this end, let us study the second 
term of the series in Eq. flT7|) . It corresponds to a one- loop diagram and is actually given by 

2 

T (2) bi] = ^Jd d xJ d d yG{y-x-m 2 )G{x-y-m 2 ) V \{x)rt{y), (22) 

which by upon substitution of the free two-point Schwinger function associated with the (fi2(x) 
heavy field and the introduction of I(p 2 ,m 2 ) as 

HP 2 , m») = 7 ± n [ d d q ^ 9 ■ ^fl" ^ , -x , (23) 

can be written as 



(2n) d J H (q 2 + m 2 )((p + q) 2 + m 2 .)' 



a 2 



T (2) bi] = J d d x J d d ytf(x)tf(y) J d d pe-^y-^m d -H{p\ml). (24) 

In the regularization and renormalization procedure we have to eliminate the poles and their 
residues adding counterterms in the Lagrangian density, consequently let us study the I(p 2 ,m 2 ). 
Using the Feynman parametrization |p9[| , it is possible to write I(p 2 ,m 2 ) as 



I{p\ ml) = N d {-- + 0(e))(l - I) I dt(^t(l -t) + (25) 



e 2 Jo m 2 

where Nd is the area of the sphere S d -i/{2n) d . The expression given by Eq. ( p5]) contains a power 
of a binomial in a flat Euclidean d-dimensional space. When d is even, the power is an integer and 
the simple use of Newton's binomial theorem will give us a very direct way of evaluating I(p 2 , m 2 ). 
When d is odd, the expansion of (1 + ■j^st(l — t))~ 2 yields an infinite power series. Since we are 
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using dimensional regularization we have an infinite power series. Note that the generalization of 
the binomial series is valid for any complex exponent p. In other words we have an everywhere 



convergent power series in p, hence a continuous function on p in the complex plane j3(J. If we 
define 

C° = 1, C 1 = — , C 2 = 21 ~ • • • (26) 



until 



r(r-l)..(r -k + 1) 



C? = -± ^ ^, (27) 

where r = | — 2, it is possible to write /(p 2 , m 2 ,) as 

J 1 oo 2fc »l 

/(p 2 , m 2 ) = (1 - + 0(e)) £ CL 2 ^ / dt(t(l - t)) fc . (28) 

Let us use the definition of Euler's integral of first kind B(a,j3) given by [j3l[ . 



B(a,P)= f dxx a -\l - xf~ l = r „,„} , Rea>0,Re/3>0. (29) 
Jo T(a + p) 

Substituting Eq. (|28|) and Eq. fl29|) in Eq.(|4]) we find that the second term of the series that 
represents the effective action, T^ 2 ) [(pi) can be written as 

r (2) M = J d d x J d d yri(x)rf(y) J d d q e^"^" 4 

j 1 oo 2A? 

((1 " M— + 0(e)) J2 CL 2 ^B(k + 1, * + 1). (30) 

/ e k=0 2 m 2 

If we choose the self-interacting part of the field <pi(x) to be Xiipf(x), it is possible to define the 
renormalized coupling constant \r subtracting the polar part. Consequently, the effective action 
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for the <fi(x) field is given by 



r e //bi] = 1 / d d x ^i(x)(-A + m 2 1 )^ 1 (x) + X R f d d x <p\{x) 



2 A 

+ TTTtV^ / d d x <pl{x){-A + m 2 M(x) + O(^) 2 . (31) 



4!(47r) 2 mi 

Note that the terms k = 3,4.. are not divergent (in a four dimensional theory) and although 
they contribute to the effective action, in the limit m 2 — > oo, the heavy field <f2(x) decouples 
from the light field y?i(x). The effect of the heavy field appears only modifying the values of 
the renormalized mass mm and the coupling constant \r. Thus we showed that the Euclidean 
version of the Appelquist-Carazzone decoupling theorem works in this specific model. Another well 
known example where the decoupling theorem can be used is in quantum electrodynamics, where 
for energies much lower than the electron mass it is possible to construct a derivative expansion of 
the Maxwell field effective action integrating out the Dirac field. This is an expected result, since 
we known that the theorem is valid for renormalizable theories without spontaneous symmmetry 
breaking or chiral fermions. The above discussion justifies the approach used by some authors that 
have been using the Euler-Heisenberg Lagrangian density to investigate the radiative correction to 
the Casimir effect [jit]] , although these radiative corrections are of no phenomenological significance 



as was pointed out by Melnikov |]14| . For a careful discussion of effective Lagrangians in quantum 
electrodynamics, see for example ref. ||32|| . 

It is important to stress that instead of obtain also a effective theory for the light field as have 
been discussed by many authors in finite temperature field theory |33j , assuming some ultraviolet 



cut-off A and integrating over the Fourier modes with wavelenght smaller than A 1 , we are assum- 
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ing that the functional integral must be taken over the space of the functions that vanish on the 
boundaries. One way to implement this is to introduce delta functions in the functional integral. 
This is equivalent to evaluate the functional integral over a space of functions that satisfy the 
boundary conditions. This is the procedure that we are addopting. It is clear that this procedure 
will introduce additional surface divergences that can be eliminated by surface counterterms, and 
in the end we have the effective model for the light modes that satisfies boundary conditions over 
some surfaces. 



3 Finite size effects and the two and four-point Schwinger 
functions in the one-loop approximation 

In the last section, we studied a very simple model of two Euclidean massive scalar fields where 
the decoupling theorem can be used after the construction of an effective theory of the light field. 
We have shown that in our model of two massive self-interacting scalar fields, the heavy modes 
associated with the ^{x) field completely decouple from the light ones associated with the light 
field (pi(x) in the limit m 2 — > oo. In the case of a \\(p\{x) self-coupling, the only effect of the 
former is a modification of the mass mi, and the coupling constant of the light field. We are 
reducing the problem in this manner since we are able to concentrate in such a one field theory, 
i.e., we will consider a \\ip A (x) self-interacting model. We will consider that the field (p(x) depends 
on d— 1 unbounded coordinates that we call f, and one bounded coordinate to which we will refer 
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to as z that will be assumed to lie in the interval [0, L\. If we exclude the possibility of periodic 
or anti-periodic boundary conditions, this choice obviously breaks the full translational invariance 
because we have to assume boundary conditions on the hyperplanes z = and z = L. 

To write the full renormalized action for the theory with boundaries we need two regulators, 
the first one being the usual e that is introduced in the dimensional regularization procedure and 
the second one that we call t] representing the distance to a boundary. According to this the full 
renormalized action must be given by |23|: 



+ 1 d d - 1 r(c 1 ( V )v 2 (r,0) + c 2 (r ] )<p 2 (r,L)) 

+ J d^ricM^ir, 0) + c 4 ( V )<p\r, L), (32) 

where A(e), B(e) and C(e) are the usual coefficients for the bulk counterterms and the coefficients 
Ci(rj) i = 1, ..4, which depend on the boundary conditions for the field, are the coefficients for the 
surface counterterms. As usual all of these coefficients must be calculed order by order in pertur- 
bation theory. We are considering two different possibilities for the boundary conditions, namely 
Dirichlet-Dirichlet {DD) and Neumann- Neumann (NN) boundary conditions. These boundary 
conditions are given respectively by 

(p(r,z)\ z=0 = (p(f,z)\z=L = 0, (33) 

and 

—^r, z )\ z=0 = —^ z )\ z=L = 0. (34) 
16 



The system we are interested is invariant only under translations along the direction parallel to the 
plates, implying that what is conserved is not the full momentum but the (d—1) dimensional paral- 
lel momentum p. For such conditions, a more convenient representation for the n-point Schwinger 
functions is a mixed (p, z) one. A Euclidean scalar field <p(x) satisfying certain homogeneous 
boundary conditions on z = and L can be expanded in Fourier series as: 

*) = 7TW E M « W / d d -'p UP) e ip ", (35) 
(2ir) 2 n J 

where p is the continuum parallel momentum, and u n (z) stands for the eigenf unctions of the 
operator —j^: 

- jgi u ^ z ) = k l u n( z ), (36) 

where k n — n — 1, 2.. for DD b.c and n = 0, 1,2.. for iVTV boundary conditions respectively. 
The main difference between both boundary conditions being the presence of the zero mode. The 
free two-point Schwinger function for the theory can be expressed in the following form: 

4V, * ?, = ^KM / ^V (p + kl + m ,y (37) 

Note that in this section we have changed the notation as follows: m\ — > m and also G(x, x'; mi) — > 
Go(x, x'). It is useful to define also p — f — r'. When considering DD boundary conditions, one 
finds that the free two-point Schwinger function is explicitly given by 

Go\p,z,z') = — £sin(^)sin(^^)/ n (L,m,rf,p), (38) 
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where 

1 r e^-P 
^ m, d, p) = j-^ J d d -i P { _ p + {If)2 + m2y (39) 

It is clear that the family of I n (L, m, d, p) functions can be thought of as the free propagators of a 
tower of massive scalar fields in (d— 1) dimensions, the effective mass of each mode being given by 
M\ = m 2 + (^) 2 . This is to be expected since our theory has been formulated in a compactified 
space. From an even simpler point of view, I n (L,m,d, p) is nothing but the Fourier transform of 
a "spherically" (SO(d — 1)) symmetric function of the parallel momentum p. 

We begin the study of the interacting theory by building the one- loop correction (G^ (Xi,x, x')) 
to the bare two-point Schwinger function G^\x,x'), for both the DD and NN boundary condi- 
tions. Using the Feynman rules we have 

G?\\ 1 ,f 1 ,z 1 ,f 2 ,z 2 ) = ^ J d d - l r J^dz G { *\r 1 -r,z 1 ,z)G ( *\d,z)Gi?\r-r 2 ,z,z 2 ). (40) 

Here we would like to point out that even though the functions G^\fi — f 2 , zi, z 2 ) and G^\f 2 — 
r 3 , z 2 , z 3 ) are singular at coincident points (fi = f 2 , Z\ = z 2 ) and (f 2 = r 3 , z 2 = z 3 ), the singularities 
are integrable for points outside the plates. Using the notation Gq'(0,z) = TDD{L,m,d, z), 
a straightforward substitution yields the order Ai correction to the bare two-point Schwinger 
function in the one-loop approximation, for the case of Dirichlet boundary conditions: 

„( 2 ) M _ v 2Ai f L ~ . nixz 1 . n-KZ . n'nz . n'irz 2 

Gi ; (Ai,n -r 2 ,z 1 ,z 2 ) = , w t 2 / dz ^ sm(— — ) sm(— ) sm(— ) sm(— — ) 

v / ^ n,n'=l 

/ e ip(f!-r 2 ) 
^ V + (f) 2 + + m 2 + ^) Tdd{1, m ' d ' zl (41) 
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where, since we are using dimensional regularization techniques, we have introduced a dimensional 

parameter /i, defining a dimensionless coupling constant A = Ai/i 4_d , and the expression for the 

amputated one- loop two-point function Tdd(L, m, d, z) is given by 

2 00 nixz r 1 
Tdd(L, m, d, z) = w ^ i ^ ( — ) j^ P [p2 + (¥)2 + ^ (42) 

In the case of Neumann- Neumann boundary conditions the expression for the amputated one-loop 
two-point function can also be found following the same procedure, and it is given by 

+^g~ 3 (™>/^ V + re)' +m ') - (43) 

Both T DD (L,m,d, z) and T NN (L,m,d, z) diverge in their continuum momenta integrals and also 
in the discrete mode summation. Using the Feynman rules, G 2 : \\, x 2 , %3, X4), i-e., the 0(A 2 ) 
correction to the bare one-loop four-point Schwinger functions, is given by 

G^\\,f 1 ,z 1 ,f 2 ,z 2 ,f 3 ,z 3 ,f A ,z 4 ) = y Jd^ 1 rJd d - 1 r f l L dzj L dzfG?\r 1 -f,z 1 ,z) 

G^(f 2 -r,z 2 ,z)(G^(r-r',z,z')) 2 
G$\? - f 3 , z', z 3 )G$\? - f 4 , z', z 4 ). (44) 

Again, all Go's are singular at coincident points, but the singularities are integrable for points 
outside the plates, except for Gq (f— r', z, z'). 

In the next section we will begin the renormalization program for the massless one-loop two- 
point Schwinger functions for the case of Dirichlet-Dirichlet boundary condition. The study of the 
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complementary set of boundary conditions, namely NN boundary conditions can be performed 
along the same lines. When the fields satisfy NN boundary conditions that infrared divergences 
for massless fields appear, and in fact, such divergences come from the zero mode contribution, 
so the two-point Schwinger function for the case of Dirichlet-Dirichlet boundary conditions is IR 
finite for m = 0. For the case of NN we must have a finite Euclidean volume to regularize the 
theory in the infrared. Another way to deal with the infrared divergences in the case of Neumann- 
Neumann boundary condition in to perform a resummation of the daisy diagrams [M\ |35[] |3B 



37| . Although this procedure is standard in the study of scalar models where the translational 
invariance is maintained, for systems where the translational invariance is broken, the problem of 
how to carry out the resummation program still remains open. 

4 The regularized one-loop two and four-point Schwinger 
functions 

In this section we would like to discuss in detail how to implement the one-loop renormalization 
program in finite size systems where flat surfaces break the translational invariance. Thus, the 
aim of this section is first to analyze the structure of the divergences associated with the one-loop 
two and also four-point function for the case of Dirichlet-Dirichlet boundary boundary conditions. 

The amputated one-loop two-point Schwinger function TDD(L,m,d, z) can be decomposed in 
a translational invariance part and another one that breaks the translational invariance, indeed 
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using algebraic identities |3T| |3§] one gets 

T DD (L,m,d,z) = fi(L,m, d) - f 2 (L,m,d,z), (45) 

where the functions fi(L, m, d) and f2(L, m, d, z) are given respectivelly by 

1 00 r 1 

h{ L,m,d) = — 

and 



x/r 1 \ 1 /"jd-i 1 cosh((L - 2z)^p 2 + m 2 ) . . 

f 2 (L,m,d,z) = — — -3— / gT 1 p^=== /_ =■ — -. (47) 

J V ' ' ' 7 2(2ix) d - 1 J y ^p 2 + m 2 senh{L^p 2 + m 2 ) V ; 

The amputated one-loop two-point Schwinger function for the Neumann-Neumann boundary con- 
ditions, T NN (L, m, d, z) can be similarly split up as 

T NN (L, m, d, z) = fx(L, m, d) + f 2 (L, m, d, z) (48) 

The above decompositions of Tdd(L, m, d, z) and TNN{L,m,d, z) have the same functional form 
and as we stated before, some of the divergences come purely from the bulk while others depend 
on the distance to the boundaries. Indeed, since fi(L,m,d) does not depend on z, it only carries 
information about the divergences on the bulk. These divergences occur not only in the discrete 
sums but also in the momentum integrations. After the identification: (3 = 2L fi(L,m,d) is for- 
mally proportional to the amputated one-loop two-point function of the theory assuming that the 
system is in thermal equilibrium with a reservoir at temperature To deal with the divergences 
of fi(L,m,d), or equivalently, the one-loop two-point Schwinger functions at finite temperature 
we have to do frequency sums and (d — 1) dimensional integrals for the continuum momenta. One 
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way to perform the integrals with Matsubara sums is to analytic extend away from the discrete 
complex energies down to real axis with the replacement of the energy sums by contour integrals 



39 1 [40]. Another way is to use dimensional regularization and afterwards to analytically extend 



the modified Epstein zeta function which appears after dimensional regularization. Direct use 
of dimensional regularization identities and the analytic extension of the modified Epstein zeta 
function in the sum given by Eq.(^) which defines f\(L, m, d), give us a polar part (size indepen- 
dent) plus a size-dependent analytic correction. The mass counterterm (the principal part of the 
Laurent series of the analytic regularized quantity) generated by fi(L,m,d) is size independent, 
because the finite temperature field theory has no temperature dependent counterterm. Observe 
that the non translational invariant part of the amputated one-loop two-point Schwinger function 
expressed by Tdd(L, m, d, z) and T^n(L, m, d, z) have the same z dependent part in modulus but 
with opposite signs. 

Since we have shown that the Tdd(L, m, d, z) and Tnn(L, m, d, z) can be split into two functions 
fi(L,m,d) and f2(L,m,d, z) and since as we have just discussed the behavior of fi(L,m,d), we 
can now turn our attention to the study of the divergences contained in f?(L, m, d, z). We begin 

d 

by an angular integration (d d ~ 1 p = p d ~ 2 dpdQd-i and, / dVtd = ) that leads to an alternate 
expression for the non translational invariant part f2(L,m,d, z), namely 

, ., , s 1 , . f°° , p d ~ 2 cosh((L - 2z)Jp 2 + m 2 ) , . 

f 2 (L, m, d, z) = -h(d) / dp Z—^ U j^JL^— L. (49) 

2 Jo yp + m z stnh{Lyp z + m A ) 

Using the change of variables s = \Jp 2 + m 2 in the above expression yields the following formula 
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for f2(L, m, d, z): 

| POO d _g 

f 2 (L,m,d,z) = -h(d) ds(s 2 - m 2 ) — C osh((L - 2z)s)(smhL S y 1 , (50) 

2 Jm 

where h(d) is an analytic function of d given by h(d) = r ^d-iy We now start studying the 

massless case following Fosco and Svaiter [ In fact, we are particularly interested in examining 
the limits (z — > + and z — > L~) which do obviously contain the information about the effects of 
the boundaries. In order to fulfill this goal we introduce two new variables x = Ls and q = zs in 
terms of which we can write f 2 (L, m, d, z)\ m=Q as 



f 2 (L,m,d,z)\ m=0 = 2 L d-2 J dxxd 3 i cotilx ~ 1) cosh(-^) 

/ dqq d - 3 e- 2q . (51) 

JO 



The second term of Eq. (^T|) give us the well known result that for a massless scalar field in d = 4 
the one-loop vacuum fluctuations diverges as \ if we approach the boundary at z = f4"2"f. The 



other term of Eq.(^Tj) should behave as ^ L _^ d _ 2 . To see this let us investigate the behavior of the 
first integral of f 2 (L, m, d, z)| m =o near the boundary at z — L. In order to do this, we make use 
of two formulas involving the definition for the Gamma function, and also another well known 
integral representation for the product of the Gamma function times the Hurwitz zeta function 
given by 

poo R 

/ dx z^e-^fcothac - 1) = 2 1 ~'T(//)C(//, - + 1) Re(/3) > 0, ReU) > 1, (52) 
Jo 2 
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where ((z, a) is the Hurwitz zeta function defined by |3~T|j 

OO J 

C(^a) = E 7 ~ w ite(«)>l, a ^0,-1, -2... (53) 

From the definition of the Gamma function and Eq.(|52~D in Eq.(|5"l"D we may write the following 
closed expression 

h(d) 



f 2 (L,m,d,z) 



m=0 



2L d - 2 
1 

(2z) 



2^V(d -2)[((d-2,- + l) + C(d - 2, -- + V 
77 ^—h(d)r(d-2). (54) 



From this last expression and using the definition of the Hurwitz zeta function giving by Eq. fl53|) 
it is evident that the regularized f 2 (L,m,d, z)\ m=0 has two poles of order (d — 2), one at z = 
and another at z = L. 

To study the massive case, from the expression given by Eq . (|50D it is possible to write f 2 (L, m, d, z) 
in a more convenient way by: 

f 2 (L,m,d,z) = f 2 i(L,m,d,z) + f 22 (L,m,d, z). (55) 

where f 2 i(L,m,d,z) and f 22 (L,m,d,z) are 

f 21 (L,m,d,z) = -h(d) ds(s 2 - m 2 )^ e~ 2z \ (56) 

2 Jm 



and 



f 22 (L,m, d, z) = -h(d) / ds(s 2 — m 2 )~(coth Ls — 1) cosh 2zs. (57) 

2 Jm 
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Using an integral representation of the Bessel function of third kind or Macdonald's functions it 
is possible to find the following closed expression of fi 2 (L, m, d, z) given by 

f 21 (L,m,d,z) = ^ 2 J. )d _ 1 (^)^K^ 1 (2mz). (58) 

For small z and finite m we have the asymptotic formula K u (z) ss 2 v ~ l Y{v)z~ v , thus for z — > + , 
the function f2i(L, m, d, z) diverges as -3=5. To calculate f 22 (L,m,d, z) we will use the same 
method that we used in section II. Again note that it contains a power of a binomial. Making use 

of the generalized binomial formula gives 

2 00 

(l-^)^ = E(-l)^(-r ? (59) 



f 22 (L,m,d,z) = f^{-l) k CU(Lmf k f° duu d -^ 2 \cothu - 1) cosh(^) (60) 

2L a z , n 2 J Lm L 



and introducing a new variable u = Ls we obtain 

fc=0 

Our next step is to show that this result can be expressed in terms of the Hurwitz zeta function. 
A natural way to achieve the proof is to split f 22 (L, m, z, d) as a sum of two terms 



f 22 (L, m, d, z) = f£(L, m, z, d) + f 22 (L, m, z, d), (61) 

where 

fc< d-3 

1 . 2 noo 111 7 

f<(L,m,z,d) = -—— - Yl C (1 \d,k){Lm) 2k duu d ^- 2k {cothu - 1) cosh(-^), (62) 

4.L J Lm L 



k=0 

and 



1 00 r°° 2n7 

f>(L,m,z,d) = --—, ]T C {1 \d,k){Lm) 2k duu d - 3 - 2k {cothu - 1) cosh(-=-). (63) 

4±/ , , J Lm Li 
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Here we have introduced C«(d, k) = (-l) k CUh(d) and also C^(d, k) = T( 2 d d Ztl k) fc) it is 

2 

possible to write Eq.(^) in the following way: 

k < d ~ z 

f<(L,m,z,d) = -— ^ c( 2 )(rf,A;)(Lm) 2fc (c(rf-2-2fc,-| + l) + C(rf-2-2A;,| + l)) 

— . V C {1 \d,k)(Lm) 2k duu d - 3 - 2fc (cothw-l)cosh(-H, (64) 



4:L d ~ , „ 

fe=0 

where the singularities of f^(L, m, z, d) appear at z — > L. Turning our attention to f^(L, m, z, d), 
it is clear that in the expression above we see that the surface divergences are the same as we 
studied before in the massless case. 

We now turn our attention back to the four-point Schwinger function in the one-loop approx- 
imation. Introducing new variables as u± = z ± z f , the two-point Schwinger function in the 
tree-level can be split into 

G$\p,z,z') = G { *\p,u + )+G®(p,u-), (65) 
where making use of the definition of I n (L, m, d, p) given by Eq.(|39|) we have 



G±\p,u±) = =p— ]T cos(^^)I n (L,m,d, p). (66) 

L n=l L 

Before continue, let us present a explicit formula of the free two-point Schwinger function in terms 
of Bessel functions. Defining an analytic function g(d) by 



i r(^ ) 



9(d) = — N ^ / 3 , (67) 
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it is possible to show that we can write G± (p, u±) as 

G<£\p, u ± ) = T 4y- £ cos{— f + m 2 ) — K^(p(m 2 + (— ) 2 ) 2 ) (68) 

P 2 £ n=l L L 2 I, 

Using Eq.(|65|) and the above formula gives us the explicit expression for the two-point Schwinger 
function in a generic d dimensional Euclidean space confined between two flat paralel hyperplanes 
where we assume Dirichlet-Dirichlet boundary conditions. It is hard to use the above expressions 
for G±\p, u±) to investigate the analytic structure of the four point function G^\\, fx, zi,r 2 , z 2 , r 3 , z 3 , r 4 , 24), 
given by Eq. (pH|) , for both the bulk and near the boundaries, Nevertheless from Eqs. (|58"D and 
Eq.([39|) it is clear that the divergences of the four-point function in the one- loop approximation 
appear at coincident points and therefore the singular behavior is encoded in the polar part of 
M(A, L, m, d) given by 

M(\,L,m,d) = A 2 J d d - l r J d^V J dz J dz' F^f 1 , z, z')^ (f - ?, z, z')) 2 , (69) 

where F(f, f", z, z') is a regular function. As with the one-loop two point function, it is not difficult 
to realize that the above equation has two kinds of singularities, those comming from the bulk 
and those arising from the behavior near the surface. As before the behavior in the bulk is as that 
found in thermal field theory and consequently we will only discuss the singularities that arise 
from the boundaries. This can be done by studying the polar part of M(A, L, m, d) given by 

M(A, L, m, d) = — [ L dz [ L dz'F{z, ^{Gf (0, z, z')) 2 , (70) 
2 jo J 
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where T{z,z'^) is a regular function. Now, we recall that the form of G± (p, u±)\ p= q is given by, 

from where it is not difficult to show that the free correlation function is given by 

G ( *\p, z, z') | p=0 = f 2 (L, m, d, ^) - f 2 (L, m, d, ^±). (72) 

For the sake of simplicity we will discuss only the massless case since the singularities of the massive 
case have the same structure as the massless one. The function f 2 (L,m,d, ^) is non singular in 
the bulk, i.e. in the interior of the interval [0, L], while f 2 (L, m, d, ^f-) has a singularity along the 
line z = z'. Indeed, closer inspection shows that for < z, z' < L the only singularities are those 
at u + = 0, u + = 2L and also U- = 0. The former two are genuinely boundary singularities (the 
two conditions imply z, z' — > or z, z' — > L) while the other comming from z = z' in the whole 
domain is just the standard bulk singularity. In fact, using the structure of the two point function 
and showing just those terms from which singularities might arise, one finds that the counterterms 
for M are given by 

" pole f dz So dz ' [ (z^j^ + (2L - z- zy=i + V^Y~> + - ]2 - (73) 

where C iy i = 1, ..3 are regular functions that do not depend on z or z'. From this discussion 
it is clear that in order to render the field theory finite, we must introduce surface terms in the 
action. This is a general statement. For any fields that satisfy boundary condition that breaks the 
translational invariance, in addition to the usual bulk counterterms, it is sufficient to introduce 
surface counterterms in the action to render the theory finite. 
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5 Boundary effects and renormalization 

In the last section we present the one-loop renormalization of the Xip 4 model, and we considered 
that the field <p(x) depends on d — 1 unbounded coordinates that we call r and one bounded 
coordinate defined in the interval [0,L]. The boundary conditions on the hyperplanes z = and 
z = L are the usual Dirichlet-Dirichlet and also Neumann-Neumann boundary conditions. 

In this section we would like to discuss briefly the global approach, used to define the Casimir 
energy associated with any field in the presence of surfaces with well defined geometric shape. The 
crucial conceptual question is the meaning of the renormalized vacuum energy associated with any 
field in the presence of any macroscopic structure that divides the space into the internal and the 
external region. It is important to keep separate different situations. In the case of the parallel 
plates, the region outside the plates is the union of two simple connected domains and both have 
the same geometry of the internal region. In this situation the Casimir renormalization procedure 
is well defined and the renormalized vacuum energy is unambiguous defined. In the case of the 
spherical an the cylindrical shell, the contribution of the exterior modes are not cancelled out in the 
Casimir renormalization procedure. It is not difficult to understand the origin of the problem, as 
has been extensively discussed in the literature. If we are assuming perfectly reflecting boundaries, 
by the Weyl theorem we know that the assymptoptic distribution of eigenvalues of some elliptic 
differential operator is related with the geometric invariants associated with the surface where 



the field satisfy some boundary condition [43] [4|]. Consequently in the regularized energy we 



have divergent terms proportional to the volume, area, etc. It the Casimir definition of the 
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renormalized vacuum energy it is not possible to cancel the area contribution for a generic surface. 
The generalization of the Weyl's expansion can be done investigating the trace of the heat kernel 
on a specified manifold with boundary. We conclude that the assumption of perfect conducting 
static boundaries with a generic shape introduces new problems in order to define the renormalized 



vacuum energy of a quantum system in the presence of these macroscopic objects |46] [47 



48[ . If someone insists in the assumption of perfect conducting boundaries there are different 
ways to solve the problem of infinite energy associated with the configuration. One is to introduce 
counterterms concentrated on the boundaries, as has been discussed by the authors that use the 



generalized zeta function method p9[ . A different approach is to smooth out the plate surface by 



a classical potential |K| Hit- It is clear that the introduction of a classical potential V(x) does 



not solve the problem of surface counterterms since in this situation we have to renormalize the 
potential. A very simple situation is the case of a background field where to compute the effective 
action we have to evaluate the the following Fredholm determinant where we are assuming that 
the positive potential is a large quantity. 

D(V) = det (-A + m 2 + 7(a;))(-A + m 2 )" 1 . (74) 

For sufficient regular but large V(x) it is possible to show that for d = 4a counterterm quadratic 



in V is required in order to eliminate the divergences of the model [52 j. Thus the introduction 
of a classical potential V(x) trying to improve the unphysical boundary condition does not solve 
the problem of surface counterterms since in this situation we have to renormalize the potential. 
Instead of smoothing the plates surfaces, a more fruitful approach to avoid surface divergences, 
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discussed by Kennedy et al j|9| is to treat the boundary as a quantum mechanical object. This 



approach was developed recently by Ford and Svaiter [53| to produce finite values for the renor- 
malized < <p 2 {x) > and other quantities that diverge as one approaches the classical boundary. 
We would like to stress that there will not be any surface divergences in a more exact treatment, 
however one can still make the case that surface counterterms are a useful phenomenological ap- 
proach for dealing with the apparent surface divergences withouth going into the complexity of 
the more exact approach. 



6 Conclusions 

In this paper we discussed the approach of effective theory to perform calculations in field theory 
in the presence of macroscopic structures. We first assumed the theory of two interacting massive 
scalar fields <pi(x) and <f2{x) with masses m\ and m 2 satisfying the condition m 2 >> mj. Inte- 
grating out the modes of the field y? 2 (x) we obtained an effective Lagrangian density for <pi(x). In 
the limit m 2 — > oo the field <p2{x) decouples from <p\{x), the only effect of (f2{x) being modifying 
both the value of the renormalized mass mi and the coupling constant of the light field (pi. Thus 
we considered the ^ipf model on a ci-dimensional Euclidean space, where all but one of the coor- 
dinates are unbounded. Translation invariance along the bounded coordinate, z, which lies in the 
interval [0, L], is broken because of the boundary conditions (BC's) chosen for the hyperplanes 
z = and z = L. Two different possibilities for these BC's boundary conditions are considered: 
DD and NN, where D denotes Dirichlet and N Neumann, respectively. The renormalization 
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procedure up to one-loop order was implemented. The main result of our investigations is that 
in the presence of boundaries where the field satisfies some boundary condition, the augmented 
action with surface counterterms can deal with the surface divergences that appear in the one-loop 
Feynman diagrams. 

There are several directions for future research in field theory in the presence of surfaces, of 
which we would like to emphasize two. The first one is to implement the renormalization program 
beyond the one-loop approximation, where overlapping divergences emerge. The second one is 
related to the infrared divergences. As we discussed, one way to deal with the infrared divergences 
in the case of Neumann-Neumann boundary condition in to perform a resummation of the daisy 
diagrams, although this procedure is standard in the study of scalar models at finite temperature, 
for systems where the translational invariance is broken, it is a open problem how to perform the 
resummation program. Both subjects are under investigation by the authors. 
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